How to calculte spectra of Rabi and related models? 
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We show how to use properly the Bargmann space of entire functions in the analysis of the Rabi 
model. We are able to correct a serious error in recent papers on the topic and develop a corrected 
method of finding the spectrum applicable also to other, more general systems. 
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Several papers devoted to determination of the spec- 
trum of the Rabi model have appeared recently, see e.g. 

and references therein for a description of various 
approaches to this problem. Let us recall that the Rabi 
model describes interaction of a two-level atom with a sin- 
gle harmonic mode of the electromagnetic field. In the 
Bargmann- Fock representation this model is described by 
the following system of two differential equations 



(z-t-A) 
(z-A) 



dz 



(1) 



where E is the energy, A is the atom-field coupling con- 
stant, and 2/x is the atomic level separation. In this rep- 
resentation, two component wave function ip — (-01, ^^2) 
is an element of Hilbert space = Jf x 3-C, where !H is 
the Bergmann-Fock Hilbert space of entire functions of 
one variable 2: e C. The scalar product in 3-C is given by 



if, 9) 
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/(z)5(z)e-l^l'd(5i(z))d(3(z)). 



An entire function f{z) belongs to J{ if it has the proper 
growth at infinity, for details, see Thus, energy E 
belongs to the spectrum of the problem, if and only if for 
this value of E equations ^ have entire solution -0 = 
{''Pi} 4^2) with the proper behaviour at infinity. However, 
for the considered system all entire solutions belong to 
5f^, see for an explanation. 

Equations ([T|) have two singular regular points at z = A 
and z = —A, while infinity is an irregular singular 
point. The system has a Z2 symmetry. It is invari- 
ant with respect to the involution r : given 
by T{'ipi,i>2){z) = {i>2{-z),^i{-z)). In other words, 
if (■!/)i(z), ■!/)2(z)) is a solution of this system, then also 
(V'2(— z), 2)) is its solution. We say that a solu- 
tion ip ~ {4'i,ip2) of ^ has parity a G { — 1,+1}, if 



t{iP) = aip. 

In [l[ it was pointed out that the the above mentioned 
symmetry plays the crucial role in determination of its 
spectrum. Later, in Q, it was shown that one can de- 
termine the spectrum of the Rabi model without any ref- 
erences to its Z2 symmetry. However, the method used 
in 3] has limited application as it is based on a result 
valid only for third order recurrence relations, see 0. 

Braak's approach is very interesting because it gives a 
chance for generalisation beyond third order recurrence 
relations. However in its original formulation it contains 
a flaw which invalidates its results concerning the spec- 
trum of the standard Rabi model. The aim of our paper is 
to correct results of [l| and to present our direct method 
which can be used to study of systems more complicated 
than the Rabi model. 

We applied our method to study a generalised Rabi 
model with broken Z2 symmetry introduced in 1] . In the 
Bargmann representation it is described by the following 
equations 



[z + X) 
(z-A) 



d0i 



d02 
dz 



{E- 



Az)V'l - ^02, 



(2) 



{E + e + Xz)^2 — fJ-^i- 



For e = it coincides with the Rabi model (jlj . We denote 
by p :— {E, X, fi, e) parameters of this model. Instead 
of E we will use also the parameter x :— E + X^. Let 
us assume that x is not a non-negative integer. Then 
the spectrum of the problem coincides with zeros of one 
function w(p) defined in the following way. Let 



where 



wip;y) ~ H,{y)H'2{y) - H[{y)H2{y), (3) 



i?i(y) :=HcunC(ao;y), (4) 
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and 



H2{y) ■■= HcunC(ai; 1 - y), 



(5) 



are the confluent Heun functions, see j?!, with parameters 
ao ■■= (a,/3,7,(5, ?7), and ai := {-a,'y, (3, -6, 6 + i]); in 
terms of (x, A, fj,, e) these parameters are given by a = 
4A^, P — -X + e, and 



-f=-l-x-e, (5 = 2(l-2e)A^ 

2?7 = 1 - 2^2 ^ (1 + x){x - 4A2) + e(l + 4A2) 



(6) 



Then, W{p) :=w{p] 1/2). 

Here we just recall the basic steps of the procedure 
which was formulated for the Rabi model ([1]). It is fully 
described in Q. 

Let us introduce new variables y := z + X, (pi{y) = 
e~'^yipi{y — A), for i=l,2. It is obvious that functions 
ipiiz) are entire if and only if functions (pi{y) are entire. 
In the new variables system ((IJ reads 



dy 

{y - 2A)^ =ix - 4A2 + 2\y)ip2 - MV'i, 
dy 



(7) 



Let us assume that for a given p = (x, A,y^), system ([7]) 
has an entire solution (p{y) = {ipi{y),(p2{y)), where 



n=0 



n=0 



To proceed further we assume also that x ^ Z, since the 
existence of solutions for integer x has long been estab- 
lished, see [§]. Then we insert the expression for ip2iy) 
into the first equation (I7|), and integrate it. We obtain 



-^nip), for n G N. 



(9) 



Function $„(y) can be determined recursively. In this 
way we obtain a solution f^y), which is locally holomor- 
phic at y = 0. Thus, we have a solution ip{z) which is 
locally holomorphic at z — —A. However, we are looking 
for an entire solution. In order to find conditions which 
guarantee that V'(-z) is entire we can assume, without loss 
of generality that it has a given parity a. Thus putting 

Ipaiz) = {lpi^cr{z),1p2,criz)) WC haVC i'l,criz) = aip2,a{-z) 

for all z e C. Hence, function 



Gcriz,p) := ip2.a{-z) - aijji^^iz), 



(10) 



must vanish identically for all z G C. In [8,], Braak claims 
that 

1. if Ga{z,p) = for a fixed but arbitrary z belong- 
ing to the intersection of domains of definiteness 
of ipi^a{z) and '>jj2,a-{—z), then it vanishes identi- 



cally. That is, if for a given p there exists G 
i:»(-A,2|A|) nD(A,2|A|) such that G„{z^,p) = 0, 
then Ga{z,p) — for an arbitrary z E C Here 
D(zo,r) C C is an open disc of radius r and the 
centre at zq. 

2. X belongs to the spectrum of the problem if and 
only if there exists z^ € D{-X,2\X\) n D(A,2|A|) 
such that Ga{Zi,,p) = 0, for a certain p = (a;, A,/i). 

Both of these statements are incorrect. There exist en- 
tire functions, e.g., sin(z), which have infinite number 
of zeros but they are not zero functions. So, the first 
claim needs at least some justification. We show that in 
the considered context it is not valid. To see this it is 
enough to make a contour graph of function Ga{z,p) on 
the plane (x, z), see Fig. [TJ One can notice that there 
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FIG. 1. The zero level contours of G+{z,p) (continuous lines), 
and G-{z,p) (dashed lines) on the {x,z) plane, for A = 0.7 
and /i — 0.4. 

exist X such that Ga{z,p) = identically in z, as well as 
that there exist isolated zeros of Ga-{z,p). An example 
is presented in Fig. [TJ This also shows that the second 
claim of Braak is not correct. 

We notice that the function g{z) :— Gcr{z,p), as a func- 
tion of z, satisfies a second order linear differential equa- 
tion. In fact, it can be checked directly that 

(z^ - X'^)g" + [z(l -2x- 3A2) - A] g'+ 

[Az(l - Az) + {x- X^f - A^ - g = 0. (11) 

So, fixing values Ga{z,p) and G'^(z,p) at a given z, we 
determine Gcr{z,p) uniquely. Hence, the spectrum of the 
Rabi problem in Braak's settings is determined by two 
transcendental equation Gcr(Zi,,p) = and G'^{Zi,,p) = 0, 
not just one Ga-{zi,,p) = 0, for a certain claimed 
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in [l[ . This is illustrated in Fig. [2l Here we just remark 
that choice z^, = made in [l| was lucky because, for 
X < 30, all zeros of Gcr{0,p) belong to the spectrum. 




FIG. 2. Graphs of G+{Zi,,p) (continuous line), and G\{Zi,,p) 
(dotted lines) for Zi, = —0.43 as functions of x, for A = 0.7 
and ^ = 0.4. The middle zero is simple (hence G+ ^ 0) and 
the correspoding x is not in the spectrum. 

However, as we show later there exists just one tran- 
scendental function such that x belongs to the spectrum 
if and only if it is a zero of this single function. This is a 
desirable property in a general method for finding spectra 
of systems described in Bargmann representation. An ex- 
ample of the Rabi model with broken Z2 symmetry, given 
in the end of paper [l] shows that a direct generalisation 
of the Braak approach is not obvious. 

In the Bargmann representation many models are de- 
scribed by a system of linear differential equations with 
rational coefficients. The system depends on a spectral 
parameter E. The problem is to distinguish all values En 
of E such that for E = En, the considered system pos- 
sesses a non-zero entire solution contained in the Carte- 
sian product of an appropriate number of Bargmann 
spaces Di. 

Here, for simplicity, we consider a system described by 
just one equation of the second order 



(12) 



where p{z) and q{z) are rational functions z, which addi- 
tionally may depend rationally on some parameters. Wc 
assume that the system possesses regular singular points 
si, . . . , Sm G C Infinity may be a regular or an irreg- 
ular singular point. Let pi and ai be the exponents at 
singular point s^. We want to formulate necessary con- 
ditions for the existence of a nonzero entire solution of 
equation ((T2)) . Hence, let us assume that (p{z) is such a 
solution. Thus, at each point Si one of exponents, let us 
say Pi, is a non- negative integer, and a local holomorphic 
solution at this point has the form 

00 

^,(z) := 5^ aW(z -.,)"■+". (13) 

n=0 

The radius of convergence of this series is not smaller than 
ri := miuj^i \si — Sj\. Coefficients a^*^ can be computed 
recursively. 



To proceed we have to make the following technical 
assumption. Namely, we require that 



Ui := D{s„n)nD{s,+i,n+i) ^ 



(14) 



for z = 1, . . . , m — 1. If none of ct^ is a non- negative inte- 
ger, then, up to a multiplicative constant, series (jl3l) are 
just Taylor expansions of entire solution ip{z) at points 
Si. This is why local solutions (pi{z) and 93^+1(2) have 
to be proportional for all z G Ui where both are well de- 
fined. Hence, for a certain nonzero ai € C, the function 
Fi{z) := ipi{z) — ai(p(z) vanishes for all z e Ui, so it van- 
ishes identically. But Fi{z), as a linear combination of 
solutions of equation (IT2|) . is its solution. So, Fi{z) van- 
ishes identically if and only if i^i(z^) — 0, and Fl{zi,) = 0, 
for an arbitrary z^, ^ U^. Eliminating from equations 
Fj(z^) = and F^{z^) — 0, we obtain the following crite- 
rion. If, with the given assumptions, equation has a 
non-zero entire solution, then 



W^{zi) := det 



iPi(Zi) (pt+l{Zi) 



0, 



(15) 



for a certain Zi G Ui, and i = 1, . . . , m — 1. These (m — 1) 
conditions guarantee that ip{z) is an entire function. If 
infinity is a regular singular point, they guarantee also 
that <p(z) belongs to the Bargmann space. 

Formulation of this method with relaxed assumptions 
and its version valid for systems of differential equations 
will be published in [lo|. 

Let us consider the Z2 symmetry-broken Rabi 
model An elimination of ip2{z) from this system 

gives one second order equation for ip{z) :— ^1(2;). This 
equation has the form p2)) with 



p{z) 



X + 2eX + zi-l + 2E + 2X^) 



9W = - 



E^ + 2zeX + A(A + z(-l + zX)) + 



(16) 



Then, we introduce variables y and v{y) setting z 
X{2y - 1), and 



v{y) :=exp(-2A2y)^i(A(2y-l)). 



(17) 



By direct computation we obtain that v{y) satisfies the 
following equation 



li + l , 7+1 

y y - 1 



w = 0, 

(18) 



where 



{a - P - + aP - - 7], 
{a + f3 + + a-f + /3-f) + 5 + T], 
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and parameters oq = {a, /3, 7, 5, rf) were defined earlier, 
see (HI). Equation p8)) is the Heun confluent equation. 
Its local holomorphic solution near y = is given by 
the Heun confluent function wo(y) — HeunC(ao;y). So- 
lution holomorphic near y = 1 is given by vi (y) = 
HeunC(ai, l—y), where ai :— (— a, 7, /3, — (5, (5+?]). Hence, 
local solutions 'fz^{z) holomorphic near singular points 
z = =pA are given by 

:=exp[Az + A2]HeunC(ao;(l + z/A)/2) (19) 

and 

Lp+{z) ■= exp[Az + A^] HeunC(ai; (1 - z/\)/2). (20) 



This research has been partially supported by grant No. 
DEC-2011/02/A/ST1/00208 of National Science Centre 
of Poland. 




These local solutions are given by expansion of one entire 
function, and hence their Wronskian 



(21) 



must vanish at arbitrary . The graph of the Wronskian 
W{p) as a function of energy for the Rabi model and the 
generalised Rabi model are shown in Fig. |31 and Fig. HI 
respectively. The resulting spectrum as seen in Fig. [5] 
reproduces that of [ij, with the benefit of not including 
possible isolated zeroes of G^. 




FIG. 3. Graph of Wronskian W(p) for p := {x,\,^,,t) 
(x, 7/10, 4/10,0). 




FIG. 4. Graph of Wronskian W{p) for p ■- (E,\,fi,e) = 
(S, 7/10, 4/10, 2/10). 
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FIG. 5. Spectrum of generalised Rabi model for /x = 0.7, and 
e = 0.2. 
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